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ABSTRACT 



We present the results of a weak gravitational lensing analysis to determine whether the stellar mass or the velocity 
dispersion is more closely related to the amplitude of the lensing signal around galaxies - and hence to the projected 
distribution of dark matter. The lensing signal on scales smaller than the virial radius corresponds most closely to 
the lensing velocity dispersion in the case of a singular isothermal profile, but is on larger scales also sensitive to 
the clustering of the haloes. We select over 4000 lens galaxies at a redshift z < 0.2 with concentrated (or bulge- 
dominated) surface brightness profiles from the ~300 square degree overlap between the Red-sequence Cluster Survey 
2 (RCS2) and the data release 7 (DR7) of the Sloan Digital Sky Survey (SDSS). We consider both the spectroscopic 
velocity dispersion and a model velocity dispersion (a combination of the stellar mass, the size and the Sersic index 
of a galaxy). Comparing the model and spectroscopic velocity dispersion we find that they correlate well for galaxies 
with concentrated brightness profiles. We find that the stellar mass and the spectroscopic velocity dispersion trace the 
amplitude of the lensing signal on small scales equally well. The model velocity dispersion, however, does significantly 
worse. A possible explanation is that the halo properties that determine the small-scale lensing signal - mainly the total 
mass - also depend on the structural parameters of galaxies, such as the effective radius and Sersic index, but we lack 
data for a definitive conclusion. 

Key words, gravitational lensing - dark matter haloes 



1. Introduction 

Galaxies form and evolve in the gravitational poten- 
tials of large dark matter haloes. The physical processes 
that drive galaxy formation cause correlations between the 
properties of the galaxies and their dark matter haloes. 
Hence to gain insight into these processes, various prop- 
erties of galaxies (e.g. colour, metallicity, stellar mass, lu- 
minosity, velocity dispersion) can be observed and com- 
pared (e.g. lSmith et!d"1l2009t [Graves et al.ll2009t) . This has 
lead to the discovery of a large number of empirical scaling 
laws, such as the Faber- Jackson relation (jFaber fc Jackson! 
Il976t) . These scaling laws help us to disentangle the pro- 
cesses that govern galaxy formation, and serve as impor- 
tant constraints for the theoretical and numerical efforts in 
this field. Although much progress has been made over the 
last few decades, many details are still unclear and warrant 
further investigation. 

One key parameter in galaxy formation is thought to be 
the total mass of a galaxy. Galaxies that have more mas- 
sive dark matter haloes than others attract more baryons 
as well, consequently form more stars which results in 
larger stellar masses. The relation between the stellar mass 
and the total mass of galaxies has been studied with ob- 



ser vat ions (e.g. Mandelbaum et all 2006t van Uitert et al 



2011 



iLeauthaud et al.ll2012l: iMore et al.l |20 1 It IWake et al 



2011 ) abundance-matching techniques (e. g. lBehroozi et al 



20101: IGuo et al ll2010t iMoster et al.ll2010D. semi-analytica l 



modeling (e.g. ISomerville et al.l 2008 : IZehavi et al 
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Keres et al 
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12012 ) 



2009; 
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and the two components are indeed found to be correlated. 
Another property of galaxies that is related to the total 
mass is the velocity dispersion, the luminosity weighted dis- 
persion of the motions of stars along the linc-of-sight within 
a spectroscopic aperture. The velocity dispersion provides a 
dynamical estim ate of the central m ass, and correlates with 
the stell ar mass dTavlor et al" 20101 ) and the total mass of 
galaxies (Ivan Uitert et al. 2011 ). 

A fundamental question that is of interest in this con- 
text is which property of galaxies is most tightly correlated 
to the total mass. This is interesting, because it shows which 
property in the centre of dark matter haloes is most inti- 
mately linked to the large-scale potential, and is therefore 
least sensitive to galaxy formation processes such as galaxy 
mergers and supernova activity that introduce scatter in 
these relations. The properties of galaxies we compare in 



1 



Edo van Uitert et al.: Stellar mass versus velocity dispersion 



this work are the stellar mass and the velocity dispersion. 
Note that there are various other observables that trace 
the total mass, and could have been used instead, but most 
of them are either expected to exhibit a large amount of 
scatter (e.g. metallicity), or they are closely related to the 
stellar mass (e.g. luminosity). 

The total mass of galaxies is not directly observable, and 
can only be determined by indirect means. An excellent tool 
to do this is via weak gravitational lcnsing. In weak lensing 
the distortion of the images of faint background galaxies 
(sources) due to the gravitational potentials of intervening 
structures (lenses) is measured. From this distortion, the 
differential surface mass density of the lenses can be de- 
duced, which can be modelled to obtain the total mass. A 
major advantage of weak lcnsing is that it does not rely on 
physically associated tracers of the gravitational potential, 
making it a particular useful probe to study dark matter 
haloes of galaxies which can extend up to hundreds of kpcs, 
where such tracers are sparse. The major disadvantage of 
weak lensing is that the lensing signal of individual galaxies 
is too weak to detect as the induced distortions are typically 
10-100 times smaller than the intrinsic ellipticities of galax- 
ies. Therefore, the signal has to be averaged over hundreds 
or thousands of lenses to decrease the shape noise and yield 
a statistically significant signal. However, the average total 
mass for a certain selection of galaxies is still a very useful 
measurement, which can be compared to simulations. 

It is important to note that the lensing signal on small 
and large scales measures different properties of dark mat- 
ter haloes. On projected separations larger than a few times 
the virial radius, the lensing signal is mainly determined 
by neighbouring structures, and therefore depends on the 
clustering properties of the lenses. Within the virial radius, 
on the other hand, the lensing signal traces the dark mat- 
ter distribution of the halo that hosts the galaxy and is 
therefore directly related to the halo mass. In this work, we 
ignore the lensing signal at large scales and instead focus 
at the signal at small scales. 

This work is a weak-len sing analogy of the analysis pre- 
sented in I Wake et al.l (|2012ft . who performed a similar study 
using galaxy clustering instead of gravitational lensing. One 
of their main findings is that the spectroscopic velocity dis- 
persion is more tightly correlated to the clustering signal 
than either the stellar mass or the dynamical mass. This 
implies that the velocity dispersion better traces the prop- 
erties of the halo that determine its clustering, that is the 
halo mass or the halo age. As the small scale weak lensing 
signal measures the halo mass, it allows us to disentangle 
the possible explanations of the clustering results. 

The outline of this work is as follows. In Section [51 
we discuss the various steps of the lensing analysis: we 
start with a description of the lens selection, then pro- 
vide a brief outline of the creation of the shape measure- 
ment catalogues, and finally discuss the lensing analysis. 
The measurements are shown in Section [3l and we con- 
clude in Section 0] Throughout the pape r we assume a 
WMAP7 cosmology (jKomatsu et al.ll201lft with cr 8 = 0.8, 
Q A = 0.73, n M = 0.27, n b = 0.046 and h 70 = H /70 
km s _1 Mpc -1 with Ho the Hubble constant. All distances 
quoted are in physical (rather than comoving) units unless 
explicitly stated otherwise. 



2. Lensing analysis 

In this study we use the ^300 square degrees of over- 
lapping area betw een the Sloan Digital Sky Survey (SDSS: 
lYork et all 12000ft and the Red -sequence Cluster Survey 2 
CRCS2: lOilbank et alJl20ll . We use the SDSS to obtain 
the properties of the lenses (e.g. stellar mass, velocity dis- 
persion), information that is not available in the RCS2. 
The lensing analysis is performed on the RCS2, because it 
is ^2 magnitudes deeper than the SDSS in r' . The increase 
in depth combined with a median seeing of 0.7", which is 
a factor of two smaller than the seeing in the SDSS, re- 
sults in a source galaxy number density that is about five 
times higher, and a source rcdshift distribution that peaks 
at z~0.7. Therefore, the RCS2 enables a high-quality de- 
tection of the lensing signal, even for a moderate number 
of lens galaxies. 



2.1. Lenses 

The SDSS has imaged roughly a quarter of the en- 
tire sky, and has measured the spectra for about one mi l- 
lion galaxies (|Eisenstein et al.l 12001b IStrauss et al.l l2002h . 
The combination of spectroscopic coverage and photome- 
try in five optical bands (u,g,r,i, z) in the SDSS provides 
a wealth of galaxy information that is not available from the 
RCS2. To use this information, but also benefit from the im- 
proved lensing quality of the RCS2, we use the 300 square 
degrees overlap between the surveys for our analysis. We 
matc h the RCS2 catalogues to the DR7 (jAbazaiian et aL 
12009ft spectroscopic catalogue, to the MPA-JHU DR7E 
stellar mass catalogue and to the NYU Value Added 
Galaxy Catalogue (NYU-VAGCl^l (iBlanton et al 



200E 



lAdelman-McCarthv et~aTll2008t iPadmanabhan et alj|200 
which yields the spectroscopic redshifts, velocity disper- 
sions, and the stellar masses of 1.7 x 10 4 galaxies. From 
these galaxies we select our lenses using criteria that are 
detailed below. 

The spectroscopic fibre within which the velocity dis- 
persion is measured has a fixed size. The physical region 
where the velocity dispersion is averaged is therefore dif- 
ferent for a sample of galaxies with different sizes and red- 
shifts. To account for this, we follow 1 ezanson 
and scale the observed spectroscopic velocity dispersion to a 
fixed size of R e /8 using cr spe c = <P oc (8.0r ap /i? e ) 066 , with 
r ap =1.5" the radius of the SDSS spectroscopic fiber, R e 
the effective radius in the r-band, and a^ cc the observed 
velocity dispersion. This correction is based on the best- 
fit relat ion determined using 40 galaxies in the SAURON 
sample (jCappellari et a l. 2006). However, the spectroscopic 
velocity dispersions provided in the DR7 spectroscopic cat- 
alogues are generally noisy for late-type galaxies. To obtain 
more robust velocity dispersion estimates for these galax- 
ies, we also predict the velocity dispersion based on quan- 
tities that are better determined following iBezanson et al.l 
(pOll : 



GAL 



^mod 



0M7K v (n)R e ' 



(1) 



with M* the stellar mass, n the Sersic index and Ky(n) 
a term that includes the effects of structure on stellar dy- 

1 http://www.mpa-garching.mpg.de/SDSS/DR7/ 

2 http:/ /sdss. physics. nyu.edu/vagc/ 
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Fig. 1. Comparison of the spectroscopic velocity dispersions to the model velocity dispersions for all galaxies with SDSS 
spectroscopy. The green triangles show the average spectroscopic velocity dispersion for bins of model velocity dispersion, 
the purple diamonds show the average model velocity dispersion for bins of spectroscopic velocity dispersion. The error 
bars indicate the scatter. The blue line shows the one-to-one correspondence. Only galaxies with a spectroscopic velocity 
dispersion error smaller than 15% have been used in the comparison. The velocity dispersions correlate well at z < 0.2, 
but at z > 0.2 the range in velocity dispersion becomes too small to assess whether this is still the case. The square of 
the correlation coefficient r 2 of the galaxies in the range 1.8 < log 10 (<7 mo d/ spC c) < 2.8 km s~ J is shown in the lower right 
corner of each panel. 



namics, and can be approximated by (|Bertin et aLl^OO"^ 

73 32 

gy < W > a 10.465 + (n-0.94) 2 + °- 954 - (2) 

The equation for <7 mo d is based on the results of 
I Taylor et al.l (|2010f) , who demonstrated that the structure- 
corrected dynamical mass is linearly related to the stellar 
mass for a selection of low-redshift galaxies in the SDSS. 

The stellar mass estimates in the MPA-JHU DR7 cat- 
alogues are based on the model magnitudes. The Sersic 
index and the effective radius in Equation (jTJ) , however, 
correspond to a different flux, i.e. the Sersic model flux, 
which is the total flux of the best fit Sersic model. This 
flux is also provided in the NYU-VAGC catalogue, and dif- 
fers slightly from the model flux. To calculate cr mo d consis- 
tently, we therefore scale the stellar mass with the ratio of 
th e model flux to the Ser sic model flux. 

iBezanson et al.l (|2011l) find that the model and the ob- 
served velocity dispersion correlate very well in the range 



60 km s _1 < a < 300 km s -1 , for galaxies in the red- 
shift range 0.05 < z < 0.07, and for a few galaxies with 
redshifts 1 < z < 2.5. The SDSS spectroscopic sample ex- 
tends to z ~ 0.5, and therefore contains many more mas- 
sive galaxies. To determine whether the velocity dispersions 
correlate well in this range too, we compare the dispersions 
for the complete SDSS spectroscopic sample in Figure [TJ 
We find that the velocity dispersions agree well, though at 
z > 0.2 the range in velocity dispersion becomes too small 
to assess whether the velocity dispersions are still corre- 
lated. This is reflected by the correlation coefficient of the 
log of the velocity dispersions of galaxies in the range 1.8 
< log 10 (<7 mo d/spec) < 2.8 km s _1 which we show in the cor- 
responding panels. 

To study whether the spectroscopic velocity dispersion 
and the model velocity dispersion agree equally well for 
different galaxy types, we split the galaxies based on their 
frac-dev parameter from the SDSS photometric catalogues. 
This parameter is determined by simultaneously fitting 
frac-deV times the best-fitting De Vaucouleur profile plus 
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Fig. 2. Comparison of the spectroscopic velocity disper- 
sions to the model velocity dispersions for galaxies with 
low concentration brightness profiles (frac-dev < 0.5) (left) 
and with high concentration brightness profiles (frac-dev > 
0.5) (right) in the redshift range < z < 0.2. For the latter, 
the dispersions agree very well, but for the former, we find 
that the spectroscopic velocity dispersion is roughly 0.1 dex 
higher than the model velocity dispersion. 



(1-frac-deV) times the best-fitting exponential profile to 
an object's brightness profile. The frac-dev parameter is 
therefore a measure of the slope (or concentration) of the 
brightness profile of a galaxy. In the following we will re- 
fer to lenses with frac-dev > 0.5 as galaxies with a sur- 
face brightness profile with a high concentration; these are 
bulge-dominated, which is typical for early-type galaxies. 
The lenses with frac-dev < 0.5 are referred to as those 
with a low concentration brightness profile; they arc disk- 
dominated as is generally the case for late-type galaxies. 
We select all galaxies with rcdshifts z < 0.2, and show the 
comparison in Figure We find that for the galaxies with 
high concentration (bulge-dominated) brightness profiles, 
the spectroscopic and model velocity dispersion agree very 
well. For those with low concentration brightness profiles, 
however, we find that the spectroscopic velocity dispersion 
is ~0.1 dex highe r than the model vel ocity dispersion. This 
is not surprising : iTavlor et abl (|2010f ) found that the rela- 
tion between the stellar mass and the structure-corrected 
dynamical mass has a weak dependence on the Sersic index, 
i.e. the ratio of the stellar mass and the dynamical mass in- 
creases with increasing Sersic index (see Figure 14 in Taylor 
et al. 2010). The offset in the relation between spectroscopic 
and model velocity dispersion for galaxies with low concen- 
tration brightness profiles is a direct consequence. It might 
be caused by the contribution of the disk velocity of spi- 
ral galaxies to the spectroscopic velocity dispersion. One 
could in principle apply a Sersic index dependent correc- 
tion, but we choose to use only galaxies with high con- 
centration brightness profiles, because there are very few 
lenses with low concentration brightness profiles in the ve- 
locity dispersion range we are interested in. As a test we 
repeated the analysis including all lenses, and found that it 
did not affect our conclusions. 

In Figure |3l we plot the spectroscopic and model ve- 
locity dispersion as a function of stellar mass. We only se- 
lect galaxies with redshifts z < 0.2; at higher redshifts, 
the range in velocity dispersions is too small to establish 
whether the correlation works well. The three lens sam- 
ples we use are indicated by the dashed lines. We select all 
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Fig. 3. Model velocity dispersion (left) and spectroscopic 
velocity dispersion (right) as a function of stellar mass. The 
dashed lines indicate the selection cuts for the lenses. 



galaxies with high concentration brightness profile with a 
stellar mass 10.8 < log(M*) < 11.5 in units of H^qMq; all 
with a model velocity dispersion 180 km s _1 < er mo d < 300 
km s _1 ; and all with a spectroscopic velocity dispersion 180 
km s _1 < fT spcc < 300 km s _1 and a relative error of < 0.15 
in fTspcc- With these criteria we select 4735, 4218 and 4317 
lenses respectively, and they form the lens samples of this 
study. 



2.2. Data reduction 

The RCS2 is a nearly 900 square degree imaging 
survey in three bands (cf , r 1 and z 1 ) carried out with the 
Canada-France-Hawaii Telescope (CFHT) using the 1 
square degree camera MegaCam. The photo metric calibra- 
tion o f the RCS2 is described in detail in iGilbank et al.l 
(pOllh . The magnitudes are calibrated using the colours of 
the stellar locus and the overlapping Two-Micron All-Sky 
Survey (2MASS), and are accurate to < 0.03 mag in each 
band compared to the SDSS. The creati on of the galaxy 
shape catalogues is described in detail in Ivan Uitert et all 
(|201 lh . We refer readers to that paper for more detail, and 
present here a short summary of the most important steps. 

We retrieve the Elixiijj processed images from the 
Canadian Astronomy Da ta Centre (CADC ) arch ivtfl We 
use the THELI pipeline (jErben et al.ll2005l 12009]) to sub- 
tract the image backgrounds, create weight maps that we 
use in the object detection phase, and to identify satellite 
and asteroid trail s. To detect the objects i n the images, wc 
use SExtractor (jBertin fc Arnoutsl I1996I) . The stars that 
are used to model the PSF variation across the image are 
selected using size-magnitude diagrams. All objects larger 
than 1.2 times the local size of the PSF are identified as 
galaxies. We measure the shapes of the galaxies with the 
KSB method ([Kaiser et alj |1995t iLuppino fc Kaiser! 11997 1 : 
IHoekstra et al.lll998lh usi ng the implementation described 
by IHoekstra et al. r [l998tl2000h . This implementation has 
been tested on simulated images as part of the Shear 
Testing Programme (STEP) 1 and 2 (the 'HH' method in 
Heymans et al. 2006 and Massey et al. 2007 respectively), 
and these tests have shown that it reliably measures the 
unconvolvcd shapes of galaxies for a variety of PSFs. 
Finally, we correct the source ellipticities for camera shear, 
an instrumental shear signal which originates from slight 
non-linearities in the camera optics. The resulting shape 



3 http:/ /www. cfht.hawaii.edu/Instruments/Elixir/ 

4 http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 
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catalogue of the RCS2 contains the ellipticities of 2.2 xlO 7 
galaxies, from which we select the subset of approximately 
lxlO 7 galaxies that coincides with the SDSS. 



2.3. Lensing measurement 

In weak lensing studies, the ellipticities of the source 
galaxies are used to measure the azimuthally averaged tan- 
gential shear around the lenses as a function of projected 
separation: 

<<*>(r) = ^, (3) 

^crit 

where AE(r) = £(< r) — S(r) is the difference between 
the mean projected surface density enclosed by r and the 
mean projected surface density at a radius r, and S cr ; t is 
the critical surface density: 



J crit 



4nG D,D u 



(4) 



with Di, D s and Di s the angular diameter distance to the 
lens, the source, and between the lens and the source, re- 
spectively. Since we lack redshifts for the background galax- 
ies, we select galaxies with 22 < m r i < 24 that have 
a reliable shape estimate (ellipticities smaller than one, 
no SExtractor flag raised) as sources. We obtain the ap- 
proximate source redshift distribution by applying identi- 
cal magnitude cuts to the photometric redshift catalogues 
of the Canada-France-Hawaii- T elescope Legacy Survey 
(CFHTLS) "Deep Survey" fields (lllbert et alJl200fift . 

To correct the signal for systematic contributions, we 
compute the shear signal around a large number of random 
points to which identical image masks have been applied, 
and subtract that from the measured source ellipticities. 
D etails on the calculation of this correction can be found 
in Ivan Uitert et al.l (|201l|) . This correction effectively re- 
moves both the impact of residual systematics in the shape 
measurement catalogues, and the impact of image masks 
on tangential shear measurements. Note that this correc- 
tion mostly affects large scales (>20 arcmin), as on small 
scales the lensing signal is generally averaged over many 
lens-source orientations causing the systematic contribu- 
tions to average out. The source galaxy overdensity near 
the lenses is found to be a few percent at most, confirm- 
ing that the lenses and sources barely overlap in redshift. 
Therefore, we do not have to correct the lensing signal for 
the contamination of physically associated galaxies in the 
source sample. 

Although neither the dark matter nor the baryonic 
component are well described by a singular isothermal 
sphere (SIS), the sum of the two components is remark- 
ably close for massive elliptical galaxies on scales smaller 
than the effective ra dius (e.g. iTreu fc Koopmansl l2004t 
iKoopmans et al.ll2009t ). The SIS signal is given by 



lt,sis(r) 



rE_ 
2r 



4rr(T len S 



DiDu 1 



D„ 2r 



(5) 



where te is the Einstein radius and oi ens the lensing 
velocity dispersion. Van Uitert et al. (2011) show that for 
galaxies with 19.5 < m r > < 21.5 the SIS profile gives an 
accurate description of the lensing signal up to ^300-500 
hjQ kpc. Based on the range of stellar masses and velocity 



dispersions of our lenses, we expect the majority of lenses 
to be central galaxies (see, e.g. van Uitert et al. 2011 
or Mandclbaum et al. 2006 for estimates of the satellite 
fraction for galaxies in these ranges) that are considerably 
more massive than the average galaxy. Hence the range 
over which the lensing signal is well described by an SIS is 
likely even larger for our lenses. 

To determine whether the stellar mass or the velocity 
dispersion is a better tracer of the amplitude of the lens- 
ing signal, we would ideally select lenses in a very narrow 
range in stellar mass, split those in a high and low velocity 
dispersion bin, and compare their lensing signals. A differ- 
ence between the lensing signal of the low and high velocity 
dispersion would indicate a residual dependence on veloc- 
ity dispersion. Similarly, we would like to select lenses in a 
very narrow range of velocity dispersion, split them in stel- 
lar mass and compare their signals. Comparing the lens- 
ing signals of these four bins would allow us to determine 
whether the stellar mass or the velocity dispersion is more 
closely related to the lensing signal on small scales - and 
hence to the projected distribution of dark matter. 

Unfortunately, we do not have a sufficient number of 
lenses for this approach. Instead, we have to select lenses 
that cover a larger range in stellar mass (and velocity dis- 
persion). We cannot simply split the lenses in velocity dis- 
persion and compare their lensing signals, because the stel- 
lar mass and velocity dispersion are correlated, and the 
high velocity dispersion bin also has a larger mean stellar 
mass. To account for this, we determine how the lensing 
signal scales with stellar mass, and remove this trend from 
the high and low velocity dispersion bins. We also deter- 
mine how the lensing signal scales with velocity dispersion, 
and remove this trend from the high and low stellar mass 
bin. If the lensing signal of galaxies strongly depends on 
the velocity dispersion, but only weakly on stellar mass, we 
expect a clear positive difference between the high and low 
velocity dispersion bins after we removed the trend with 
stellar mass. At the same time, we should see only a very 
small difference between the high and low stellar mass bin 
after removing the trend with velocity dispersion. Hence 
by studying the differences in the residual lensing signals, 
we can tell which observable is more closely related to the 
lensing signal of galaxies. 



3. Results 

To study whether the lensing signal mainly depends 
on stellar mass or velocity dispersion, we first have to 
determine how the lensing signal scales with these observ- 
ables. We discuss how this is done for the model velocity 
dispersion; for the spectroscopic velocity and the stellar 
mass, we follow a similar approach. The general procedure 
is summarized below. 



• We sort the lenses in model velocity dispersion, and di- 
vide them in five quintiles; 

• We measure the lensing signal of each quintile, to which 
we fit an SIS profile on scales between 50 ft^kpc and 
1 /zfoMpc. This is roughly the range where the galaxy 
dark matter halo dominates the lensing signal. This re- 
sults in five best-fit lensing velocity dispersions, aiens! 
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Table 1. Power law parameters 



300 




Fig. 4. Best-fit lensing velocity dispersion as a function of 
spectroscopic velocity dispersion (left panel, black), model 
velocity dispersion (left panel, orange) and stellar mass 
(right panel). Dashed lines indicate the best- fit linear rela- 
tion between the observable and cri en s- The linear relations 
are used to remove the dependence of the lensing signal on 
these observables. 

• We use the five values of criens to fit the linear relation 
cions = fflmod x (c mo d/200km s _1 ) + 6 mod . We show the 
measurements and the fit in Figure 01 and give the best- 
fit parameters in Table [TJ 

• We determine the median stellar mass of these lenses, 
and divide them into a low and high stellar mass sam- 
ple. We measure the lensing signal of both samples, and 
show them in the top-left panel of Figure [SJ 

• For each lens in the low and high stellar mass sample, 
we use the model velocity dispersion to calculate oi cns 
using the linear relation, and subtract their SIS profiles 
from the lensing signal. The residuals are shown in the 
middle-left panel of the same figure; 

• Finally, we determine the difference between the resid- 
ual lensing signal of the high and low stellar mass bin, 
(5(AS — AStrcnd), which is shown in the bottom-left 
panel. 

When we subtract two SIS profiles with different am- 
plitudes from each other, the result is also an SIS pro- 
file. Therefore, to quantify the residuals, we fit an SIS to 
(5(AS — AEtrond) on the same scales, and determine the 
residual Einstein radius, r^ 5 . These values can be found in 
Table [2] 

Similarly, we determine the dependence of the lens- 
ing signal on spectroscopic velocity dispersion and stel- 
lar mass. For the spectroscopic velocity dispersion, we fit 
cions = a sp oc x (cr spoc /200km s -1 ) + & spcc and for the stel- 
lar mass, we fit cti ci1s = a stc i X log(M*/10 n h^M Q ) + b ste \. 
The best- fit parameters are shown in Table Q] These trends 
are removed from the lensing signals, and the residuals are 
shown in Figure [5] (middle panel) . 



O mod 

[km/s] 


[km/s] 


^mod 

[km/s] 


180 < cr mod < 300 
100 < cr mod < 400 


44 ±48 
93 ±25 


141 ± 54 
80 ±25 


t^spec 

[km/s] 


Q^spcc 

[km/s] 


^spec 

[km/s] 


180 < (Tepee < 300 

100 < o- spoc < 400 


176 ± 43 
129 ± 21 


-5 ±50 
45 ± 23 


log(M.) 


[km/s] 


6stel 

[km/s] 


10.8 < log(Af») < 11.5 
10.5 < log(M») < 12.0 


134 ± 36 
118 ±22 


179 ±8 
178 ±6 



Notes. Best-fit power law slope a x and offset fe x that describe 
the relation between galaxy property x ('mod', 'spec' and 'stel' 
for model velocity dispersion, spectroscopic velocity dispersion 
and stellar mass, respectively) and the lensing velocity disper- 
sion for lens galaxies in the range that is indicated in the first 
column. Details of the fitting are described in the text. 



Table 2. Residual Einstein radii 



removed 
trend 


residual 
dependence 


~.res 
T E 

[/i^kpc] 


r E 

[ftrokpc] 


^mod 

M, 

<^~spcc 

M, 


M t 

^mod 

M, 

t^spec 


0.88 ±0.25 
-0.18 ± 0.24 
0.30 ± 0.25 
0.37 ±0.24 


(0.78 ±0.25) 
(-0.12 ±0.24) 
(0.42 ±0.25) 
(0.42 ±0.24) 



Notes. The residual Einstein radius, obtained by fitting an SIS 
profile to <5(AE — AEtrend) between 50 h^kpc and 1 /i^Mpc 
for a mean lens redshift of z = 0.13. The bracketed values in 
the fourth column show the results for a different linear relation 
between the galaxy property and tri ens , as detailed in the text. 



In the bottom panel of the first column of Figure [5j we 
observe that after we have removed the lensing signal de- 
pendence on model velocity dispersion, <5(A£ — AE tre nd) is 
still positive on small scales, and therefore the lensing sig- 
nal has a residual dependence on stellar mass. This residual 
dependence implies that the lensing signal still depends on 
stellar mass after removing its dependence on model ve- 
locity dispersion. In the panel next to it, where we have 
removed the dependence on stellar mass, we find that the 
difference between the residuals of the model velocity sam- 
ples is consistent with zero, i.e. the lensing signal shows no 
residual dependence with model velocity dispersion. These 
trends are reflected by the values for r^ s in Table [2] The 
third and fourth columns of Figure [5] show that if we re- 
move the dependence on spectroscopic velocity dispersion, 
the difference of the residual signal of the high and low stel- 
lar mass sample is consistent with the difference between 
the residual signal of the high and low spectroscopic veloc- 
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180<CT mod <300 10.8<log(M.)<1 1 .5 1 80<<7 spec <300 1 0.8<log(M.)< 1 1 .5 
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Fig. 5. In the top row, we show the lensing signal AS as a function of physical distance from the lens, for the lens 
samples that have been split by the median value of one of the obscrvablcs, as indicated in the plots. Red triangles 
(green diamonds) indicate the signal of the lenses with larger (smaller) stellar masses/ velocity dispersions. In the middle 
row, we show the lensing signal of the same samples after we subtracted the trend with the observable that is indicated 
on top of each column. The difference between the residual trends for the two lens samples are shown in the bottom 
row. A residual trend indicates that the lensing signal has a residual dependence on the observable indicated inside the 
corresponding panel of the first row, after removing the dependence on the observable indicated on top of that column. 
The dotted lines show the best-fit SIS profiles to the difference between the residuals. 



ity dispersion samples after we removed the dependence on 
stellar mass. Note that even though the values of the resid- 
ual Einstein radii change somewhat if we limit the analysis 
to smaller scales, which is likely mostly due to the larger 
statistical errors and because we probe different regions of 
the haloes, the trends we find do not change qualitatively. 
A quantitative characterization of the scale dependence of 
the signal could be performed with upcoming lensing sur- 
veys. 

Since we find that the lensing signal still depends on 
stellar mass after we remove its dependence on model ve- 
locity dispersion, but it does not depend on model veloc- 
ity dispersion once wc have removed the trend with stellar 
mass, our results suggest that the stellar mass is a bet- 
ter tracer of the lensing signal of galaxies than the model 
velocity dispersion. Furthermore, the stellar mass and the 
spectroscopic velocity dispersion trace the lensing signal 



equally well, as the residual Einstein radii are consistent. 
As a consistency check, we have also looked at the resid- 
ual dependence on model velocity dispersion after removing 
the trend with spectroscopic velocity dispersion, and vice 
versa. These trends confirm our previous findings: the spec- 
troscopic velocity dispersion is more sensitive to the lensing 
signal of galaxies than the model velocity dispersion. 

There is a weak indication that the lensing signal has 
a residual dependence on stellar mass after we remove 
the trend with spectroscopic velocity dispersion, and vice 
versa. This would imply that both the stellar mass and 
the velocity dispersion contain independent information on 
the projected distribution of dark matter around galaxies. 
Unfortunately, we do not have sufficient signal-to-noise to 
obtain a clear detection. 

The results depend on the linear relations we have fit 
to remove the dependence on the observables. To study 
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log(Mj[h 7D 'M ] CT mod [km/s] CT s P e C [ km / s ] 

Fig. 7. The effective radius as a function of stellar mass (left), model velocity dispersion (middle) and spectroscopic 
velocity dispersion (middle). The black dots show all galaxies with high concentration brightness profiles with z < 0.2, 
the red dots are the lenses that satisfy 180 km s _1 < cr mo d < 300 km s _1 , 180 km s _1 < er spcc < 300 km s _1 , 
^Cspec/cspcc < 0.15 and 10.8 h^M Q < log(A/„) < 11.5 /i^'Mq. By using all selection criteria simultaneously, we exclude 
large galaxies at low stellar masses and high model velocity dispersions, and small galaxies in the velocity dispersion 
samples. These selection biases could bias the lensing analysis if the lensing signal of a galaxy also depends on the effective 
radius. 



how sensitive the residual trends are on these relations, we 
have also fit them using all galaxies with high concentration 
brightness profiles in the range 100 km s _1 < a spcc < 400 
km s _1 , 100 km s _1 < cr mod < 400 km s _1 and 10.5 
/i^Mq < log(M*) < 12 ftfo M , respectively. The best- 
fit parameters of these fits are shown in Table [1] We re- 
peated the analysis using these values, and show the resid- 
ual Einstein radii between brackets in Table [5] We find 
that this docs not significantly change the results, i.e. the 
model velocity dispersion traces the lensing signal of galax- 
ies worse than either the stellar mass or the spectroscopic 
velocity dispersion. 

It is somewhat surprising that cr mo( j is a poorer tracer of 
the total mass than cr sp0 c, particularly because we observe 
in Figure Q] that they correlate well. A possible explana- 
tion is that they both trace the lensing velocity dispersion 
equally well on average, but with a different intrinsic scat- 
ter. To test this, we dr aw 10 6 galax i es from the velocity dis- 
persion function from lSheth et al.l (I2003T ) . which we adopt 
as the lensing velocity dispersion o"i ens . Next we assign them 
a spectroscopic and model velocity dispersion by drawing 
from a normal distribution with mean oi ons , using a larger 
scatter for cr mo d than for <7 spoc . Then we determine the av- 
erage cions hi the five velocity dispersion bins, as in Figure 
21 and we compare the distribution of er mo d and a spec as in 
Figure [5J Although we can choose the scatter of cr mo d and 
fj spoc such that we either reproduce Figure [5] or Figure EJ 
we cannot find a combination that reproduces both simul- 
taneously. We also tested different distributions for a mo d 
and (Tspoc instead, and found similar results. Hence a differ- 
ence in intrinsic scatter can at best only partly explain the 
different performance of cr mod and cr spcc as tracers of the 
lensing signal. 

Another option is that there is another intrinsic prop- 
erty of galaxies with some additional dependence on the 
lensing signal. To study where the samples differ, we plot 
the average effective radii and Sersic indices as a function 
o"mod and er spcc in Figure [H] We find that at low velocity 
dispersions, the values of r e and n are similar, but at high 



velocity dispersions, the lenses in the <7 spC c samples have 
larger effective radii, whilst the lenses in the cr mo d have 
larger Sersic indices. Hence the difference between the per- 
formance of <7 mod and (7 S pec could be due to an additional 
dependence of the lensing signal on the structural parame- 
ters of the lenses. 

To test whether the lensing signal depends on the size 
of galaxies, we select the lenses from the model velocity dis- 
persion sample, and remove the lensing signal dependence 
on <7 mo d- We then determine the median effective radius, 
split the lenses into a low and high effective radius sam- 
ple and measure their residual lensing signal. As before, 
we measure the difference between the residual lensing sig- 
nals of the high and low effective radius sample, to which 
we fit an SIS profile. We find r£ s = 0.53 ± 0.25 kpc, 
which suggests that the lensing signal depends on the size 
of a galaxy. However, in Figure H2 and [7] we find that the 
effective radius is correlated with stellar mass, so part of 
this residual may be caused by the dependence on stellar 
mass. Therefore, we repeat the test using the lenses from 
the stellar mass sample, and remove the lensing signal de- 
pendence on stellar mass. Then we determine the median 
effective radius, split the lenses into a low and high effec- 
tive radius sample and measure the difference between their 
residual lensing signals. We find r r ^ s = 0.14 ±0.24 kpc. 
Studying the residual dependence on Sersic index, we find 
r£f = 0.04 ±0.25 h^Q kpc for the model velocity dispersion 
sample, and r T £ s = 0.11 ± 0.24 h^Q kpc for the stellar mass 
sample. These results do not provide conclusive evidence 
that the small-scale lensing signal depends on these struc- 
tural parameters. 

Although the three lens samples overlap, they are not 
identical. Hence part of the trends we observe might actu- 
ally be due to differences in the lens samples. To test this, 
we could define a fourth sample by selecting galaxies that 
pass all selection criteria, i.e. 180 km s _1 < cr mod < 300 km 
s , 180 km s _1 < <7 spcc < 300 km s , <5cr spec /<7 spe c < 0.15 
and 10.8 h^M Q < log(M*) < 11.5 h^M Q . However, if we 
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Fig. 6. Mean effective radius and Sersic index as a func- 
tion of spectroscopic velocity dispersion (left, black filled 
diamonds), model velocity dispersion (left, black open tri- 
angles) and stellar mass (right, black diamonds). In red, 
we show the averages for the lens samples that simultane- 
ously satisfy 180 km s _1 < er mo d < 300 km s -1 , 180 km 
s _1 < (T spoc < 300 km s _1 , <5cr S pec/cspec < 0.15 and 10.8 
K^qMq < log(M*) < 11.5 h^MQ. By simultaneously ap- 
plying all selection criteria the average sizes and Sersic in- 
dices of the samples change, which shows that we implicitly 
exclude galaxies from a certain area of structural parameter 
space. 



simultaneously select on stellar mass and model velocity 
dispersion, we implicitly also select on effective radius and 
Sersic index. This is demonstrated in Figure [6] and Figure 
where we show the mean effective radius and Sersic in- 
dex for the lens samples. When we select lenses that pass 
all selection criteria, we exclude lenses with large effective 
radii and small Sersic indices at low stellar mass, lenses 
with small effective radii and small Sersic indices at low 
spectroscopic velocity dispersions, and lenses with small ef- 
fective radii and large Sersic indices at low model velocity 
dispersions. If the lensing signal of a galaxy also depends on 
its structural parameters, the lensing measurements of this 
fourth sample could be biased, making the results harder 
to interpret. 

I Wake et al.l (l2012h find that the ratio of the projected 
correlation functions of a high and low spectroscopic veloc- 
ity dispersion sample at a fixed stellar mass is of the order 
1.5-2 on scales <lh~ l Mpc, whilst the ratio of a high and 
low stellar mass sample at a fixed velocity dispersion is close 
to unity. Our results show, however, that the lensing signal 
on small scales, which mainly depends on the halo mass, 
is traced equally well by the stellar mass and the spectro- 
scopic velocity dispersion, and that potential differences are 
smaller than our statistical errors. The explanation could 
be that the difference in the amplitudes of the correlation 
functions are not mainly determined by how well they trace 



the halo mass. lWake et all (|2012h mention two other possi- 
ble causes: the relation between the velocity dispersion and 
halo age is tighter than between stellar mass and halo age, 
or tidal stripping of satellites which leads to a reduction 
in stellar mass, but does not affect the velocity dispersion. 
Although our results favour the latter two explanations, we 
cannot draw firm conclusions due to the low number of lens 
galaxies that could be used. 

It is important to note that we only study galaxies with 
high co ncentration b r ightne ss profiles, whilst the main re- 
sults of I Wake et al.l (120121) are based on a sample with 
mixed galaxy types. IWake et"aT1 (|2012f) do separate their 
sample based on colour and on morphology, and find some 
residual dependence on colour, but not on morphology, 
which suggests that their results would not have changed 
by much for a galaxy sample selected with similar criteria 
as in our work. We note, however, that the range of Sersic 
indices of our galaxies is limited, as is shown in Figure [51 
Since Sersic index and velocity dispersion generally corre- 
late well, it might be that this is diluting the effect in our 
observations, which then would support the view that the 
spectroscopic velocity dispersion is a better tracer of the 
halo mass than the stellar mass. This could be investigated 
in more detail by repeating this analysis with a larger lens 
sample. 

4. Conclusion 

In this work, we study which property of galaxies is 
most tightly correlated to the weak gravitational lensing 
signal on small scales for a sample of ^4000 galaxies with 
high concentration brightness profiles (frac-dev > 0.5) 
at z < 0.2. The properties we compare are the stellar 
mass, the spectroscopic velocity dispersion and the model 
velocity dispersion. We find that the lensing signal of 
galaxies is equally well traced by the stellar mass and 
the spectroscopic velocity dispersion. There is a weak 
indication for a residual dependence on stellar mass after 
removing the trend with spectroscopic velocity dispersion, 
and vice versa. This suggests that both tracers contain 
independent information on the projected distribution of 
dark matter around galaxies. Unfortunately, the signal-to- 
noise of our lensing measurements is not sufficient to make 
a definite statement. 

The model velocity dispersion traces the lensing signal 
significantly worse, which is surprising as the spectroscopic 
velocity dispersion and model velocity dispersion correlate 
well for our lenses. We find that this cannot be solely ex- 
plained by assuming a larger intrinsic scatter of the model 
velocity dispersions compared to the spectroscopic ones. 
At high velocity dispersions, however, the lenses in the 
Cmod-sample have smaller effective radii and larger Sersic 
indices than those in the (7 spcc -samplc. This suggests that 
these structural parameters contain additional information 
on the projected distribution of dark matter around 
galaxies. To test this, we measure how the lensing signal 
depends on the size and Sersic index of the lenses. We 
do not find conclusive evidence for a residual dependence 
on these structural parameters, which could be due to 
insufficient signal-to-noise caused by the relatively small 
lens sample of this study. 

The lensing signal on small projected separations from 
the lenses mainly depends on the halo mass. Our results 
therefore suggest that the stellar mass and spectroscopic 
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velocity dispersion trace the halo mass equally well, but 
the model velocity dispersion does worse. However, at 
larger separations, neighbouring structures contribute to 
the lensing signal as well, and we cannot exclude the 
possibility that differences between the satellite fractions 
and large-scale clustering properties of the lens samples 
also have some effect. 

Ideally, one should also remove the potential lensing 
signal dependence on the structural parameters of galaxies, 
i.e. split the lens sample both in velocity dispersion and 
structural parameters, and study the residual dependence 
on stellar mass. With the current data, we do not have suffi- 
cient signal-to-noise to perform such an analysis. However, 
we expect that at low-redshift (z < 0.2) improvement is 
possible by repeating this analysis on the complete SDSS, 
while at higher rcdshifts the overlapping area between the 
RCS2 and the data release 9 (|Ahn et al-feoll of the SDSS 
could be used. Ultimately, one could simultaneously fit all 
these parameters, i.e. A//j = f(a, M* 7 r e , n, ...), which could 
also contain products of the parameters such as crM*, and 
determine the covariancc matrix between the coefficients. 
The relative magnitude of the coefficients would give new 
insights into which observables are important, and hence 
would provide valuable insights into galaxy formation 
processes. The expec ted lensing measure ments from the 
space- mission Euclid (|Laureiis et al.|[201ll) would be ideal 
for such a study. 
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